Large scale numerical investigation of excited states in poly(phenylene) 
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A density matrix renormalisation group scheme is developed, allowing for the first time essentially 
exact numerical solutions for the important excited states of a realistic semi-empirical model for 
oligo-phenylenes. By monitoring the evolution of the energies with chain length and comparing them 
to the experimental absorption peaks of oligomers and thin films, we assign the four characteristic 
absorption peaks of phenyl-based polymers. We also determine the position and nature of the 
nonlinear optical states in this model. 



The phenyl-based conjugated polymers have at- 
tracted immense interest from physicists and chemists 
as a result of the discovery of electroluminescence in 
poly(phenylenevinylene) (PPV) [Q. In particular, a great 
deal of theoretical and experimental effort has been de- 
voted to understanding the important neutral excita- 
tions driving nonlinear optical processes, the positions of 
triplet states and the nature of exciton binding and decay 
Correlation between 7r-electrons has proved to be 
important in obtaining an accurate description of excited 
states in conjugated systems. It is thus desirable to have 
accurate numerical solutions of correlated electron mod- 
els for these systems (without recourse to uncontrolled 
approximation schemes such as configuration interaction 
approximation schemes, or perturbative schemes starting 
from a noninteracting molecular orbital or fc-space band 
description). In this paper we provide, for the first time, 
accurate results for a realistic semi-empirical model for 
poly(phcnylenc) and its oligomers, a model system for 
phenyl-based conjugated polymers. 

As our starting point for describing poly(phenylene), 
we adopt the Pariser-Parr-Pople (P-P-P) model |7j], 
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where <> represents nearest neighbors, Ci a destroys a 7T- 
electron on conjugated carbon atom i, n irj = c| Cj a and 
rii = riii + nil- We use the Ohno parameterisation for 
the Coulomb interaction, Vij = U/^Jl + (Ur^ /14.397) 2 , 
where is the inter-atomic distance in A and U = 10.06 
cV Q. The transfer integrals (tjj) and bond lengths are 
2.539 eV and 1.4 A respectively for phenyl bonds, and 
2.22 eV and 1.51 A respectively for single bonds ||. In 
the following we consider oligophenyls of N repeat units 
(phenyl rings), as depicted in Fig. ||. 



When solved accurately, the P-P-P theory has proved 
remarkably accurate in describing many of the low-lying 
excitations in a wide range of conjugated molecules ^|,^|. 
Given such debates as to the nature of exciton binding in 
PPV and the related issues of interpreting various spec- 
tral features it is intriguing as to what the P-P-P 
theory predicts for large phenyl systems. Exact diagonal- 
isation is currently restricted to systems with two phenyl 
rings. The advent of the density matrix renormalisa- 
tion group (DMRG) method |l(J has been fortuitous in 
that it has the potential to provide effectively exact nu- 
merical solutions of the P-P-P theory for systems with 
hundreds of conjugated carbon atoms. The DMRG has 
obstensively been restricted to systems with acetylene- 
like (simple, linear chain) geometry |Il]Jl2Tl to date. This 
is because it has not been possible to keep a complete 
basis for the phenyl repeat unit (4096 basis states) as 
well as sufficient system and environment block states. 
In this work we develop an accurate method for treating 
the phenyl block based on local Hilbert space optimisa- 
tionHQ. 

Our approach starts with the assumption that an op- 
timised local Hilbert space prepared in a small system 
should also be suitable for a large system. We employ 
the following procedure: 

1. The P-P-P model is solved exactly for the 12-site 
(N = 2) biphenyl molecule. The right-hand ring 
acts as the environment for the left-hand ring. A 
reduced density matrix for the left-hand ring is con- 
structed by performing a partial trace of the target 
state(s) projection operator over the states of the 
right-hand ring. The mj optimal states for the left- 
hand ring (i.e. the optimised local Hilbert space) 
are those density matrix eigenstates with the high- 
est eigenvalues. These are the natural many body 
states of the phenyl ring at the end of an oligo- 
phcnylcnc molecule, in analogy to the natural one- 
electron orbitals of quantum chemistry. 

2. For phenyl rings added to the middle of the lat- 
tice during the DMRG iterations (i.e. rings with 
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neighbours to the left and right), we apply peri- 
odic boundary conditions to the biphenyl molecule 
to again construct an optimal basis, of which the 
tub most important states are retained. We note 
that this method of determining a suitable repeat 
unit Hilbert space from a small, exactly diagonal- 
isable system is employed at the first stage of the 
"four-block" DMRG method @. 

3. Once the local Hilbert space has been optimised 
the usual inifinite lattice DMRG approach is used 
to build up oligophenyls. We use a three-block in- 
finite lattice method, the environment block being 
the reflected system block. This enables us to ap- 
ply the full symmetries of the Hamiltonian, which 
are the two planes of reflection, and particle-hole 
and spin- flip symmetries. We use the symbols A g 
and B\ u to denote states that are symmetric under 
long-axis reflection, the A g states being symmetric 
under short axis reflection, the B\ u states being an- 
tisymmetric under short axis reflection (polarised 
along the molecular axis). We use the symbol B^ u 
to denote states that are polarised perpendicular 
to the molecular axis, i.e., states which are sym- 
metric under long axis reflection and antisymmetric 
under short axis reflection. The spin-flip symme- 
try is used to distinguish singlet and triplet states, 
i.e., singlet states in the A g sector are denoted 1 A g 
whilst triplet states are denoted 3 A g . Finally, the 
particle-hole symmetry of states is denoted by a 
"+" or "— " superscript. For example, the second 
lowest state in the A g singlet sector with positive 
particle-hole symmetry is denoted 2 1 A+. 

4. Oligophenyls with an even number of phenyl rings 
are constructed by performing DMRG iterations 
with no repeat unit (i.e., just with a system and 
environment block). 

We embellish the basic approach described above with 
some further optimisations that serve to accelerate the 
convergence of energies with ttia and tub- Fig. |l| shows 
the construction of the initial phenyl bases: the end 
phenyl unit basis (for the first system block) is derived 
from the exact diagonalisation of biphenyl with open 
boundary conditions; the initial middle phenyl unit ba- 
sis is derived from exactly diagonalising biphenyl with 
periodic boundary conditions. However, before using the 
middle block basis in DMRG iterations it pays to perform 
an improvement step using the trimer (N — 3). Namely, 
we diagonalise the timer using a relatively small number 
rrijp of states in the system block, and a large number of 
states nig in the middle block. A reduced density matrix 
is then formed for the middle unit (by performing a par- 
tial trace over the system and environment block Hilbert 
spaces on an appropriate projection operator) and diag- 
onalised to give the basis used for middle units in the 
subsequent DMRG iterations. 



Fig. H shows a DMRG iterative step in the infinite 
lattice algorithm. 

1. First, a system block and its reflection are com- 
bined without a middle unit in order to treat sys- 
tems with an even number of repeat units. In this 
calculation states are retained in the blocks. 

2. Next, a middle unit is inserted in the superblock. 
However, we perform two further basis optimisa- 
tions to improve the middle block and system block 
bases in the context of this system size before ob- 
taining our final energy estimates. 

3. In the first step we improve the middle block basis 
(obtained from the trimer calculation depicted in 
Fig. nj) by again retaining a relatively small number 

(2) . 

m A states in the system and environment blocks 

(2) 

and a relatively large number m B states in the 

middle block, with < m B . A reduced density 
matrix is again formed for the middle block and 
diagonalised to give an improved middle block basis 
for this system size. 

4. A similar optimisation is then carried out for the 
system block, this time keeping a relatively small 

(3) 

number of states m B in the middle unit and a 

(3) 

relatively large number of states m A in the system 
and environment blocks. A reduced density matrix 
is formed for the system block by tracing out over 
the middle unit and environment block spaces. 

5. Finally, the energy estimates are obtained for this 

(3) 

system size by retaining tua < Tn\ states from 

(2) 

the reoptimised system block basis and mj< m B 
states from the optimised middle block basis. It is 
also at this stage that the system block and middle 
unit arc augmented to produce the system block 
used in the next iteration. 

Before discussing our results, we assess the convergence 
of the DMRG procedure with mA and ms- We first con- 
sider the noninteracting limit (U = 0) where comparisons 
can be made with exact results. Table | shows the ener- 
gies obtained for a number of states in various symmetry 
channels for N = 11 phenyl rings using tua = 280 and 
vtlb = 468. In this example — 112 and to^ = 606 
and the extra middle unit and system block optimisations 
((b) and (c) in Fig. ||) are not used. The accuracy varies 
with the state studied but even the most inaccurate state 
is resolved to within 0.02 eV. These results are a number 
of orders of magnitude more accurate than a previous ef- 
fort to apply the DMRG to the oligo-phenylene geometry 

©■ 

In the interacting case of interest the convergence of 
energies is monitored as a function of both rriA and m b ■ 
We first consider biphenyl (N = 2) in order to get an idea 
of the efficiency of the phenyl bases. In Table || we list 
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the energies of the two lowest 1 A+ states as functions of 

m } , the number of system block states retained in step 
(a) of the procedure depicted in Fig. 0. The indication 
is that even around 100 states around 3% of the total 
basis gives a reasonably accurate truncated basis for the 
phenyl block. We next consider the initial trimer step 
(N = 3) depicted in Fig. |l|. In Table [Tfi we list results 
for the 2 1 A+ transition energy as a function of m A and 

rrig\ The energy is converged to within around 0.01 
eV. We have found that, once this initial reoptimisation 
of the middle phenyl block has been done for the trimer 
system, the parameter that controls the convergence (the 
parameter to which the energies are most sensitive) is the 
usual DMRG truncation parameter tua, the number of 
states retained in the system and environment blocks. 
A typical example is given in Table IV where the 2 1 /1+ 
energy state is listed for various tua and m b ■ The value 
of is set at 140 throughout, the value of is 
varied with tub as shown, and the extra middle unit and 
system block optimisations (steps (b) and (c) from Fig. |^) 
are not used. The results indicate that the 2 1 J 4+ energy 
is resolved to within around 0.02 eV or better. 

In production calculations we use all the embellish- 
ments depicted in Fig. ||. A typical example of the pa- 

r?4 0) = 1000, 

(3) 

m A 



rameters used is given by: m ( l ] = 168, 
TO « = 340, = 140, rn§ = 740, m)£> = 290, 

m B ' = 300, rriA = 180 and tub = 520. We have not 
carried out a comprehensive study of the effect of the 
embellishments (b) and (c) but, from running a few sets 
of parameters for the same target states, and from ob- 
serving the change in energy obtained from using a re- 
duced, reoptimised basis for either the system block or 
the middle phenyl unit, the indication is that they al- 
low us to work with effectively larger values of uia and 
tub- That is, using the embellished procedure with the 
above parameter set is probably equivalent to the simpler 
scheme (without steps (b) and (c)) with uia ~ 290 and 
rns ~ 740. We are fairly confident that, in the calcula- 
tions in this paper, the lowest lying eigenstates in a given 
symmetry channel are resolved to within around 0.02 eV 
or better @. 

We now turn to a discussion of our results. We first 
embark on an investigation of the four characteristic 
absorption peaks of poly(phenylene) oligomers [ fj^]l8| 
and thin film polymers [ p^7|p9|j20[ | . Since the excitations 
of biphenyl (TV = 2) are relatively well understood, a 
knowledge of their evolution with chain length enables 
us to identify the key excitations in poly-phenyls. As 
the oligophenyls possess Z?2/i symmetry, the dipole-active 
states are polarised either along the long axis i}B^ u ) or 
short axis ( 1 i?^ u ). Fig. || shows the A^-dependence of 
transition energies of some key states, along with exper- 
imental results for biphenyl (N — 2) (see ||), oligomers 
(N = 3, . . . , 6 Jl7|arid N = 6 g|]) and polymer thin 
films (N = oo)^^,^9|j20). We note that the particle- 
hole dipole- forbidden state (l 1 ^^) lies below the dipole- 



active l x B^ u state in biphenyl. However, the l 1 ^^ state 
very weakly hybridises, and thus its energy is almost in- 
dependent of chain length, converging to 4.4 eV. This 
energy agrees very well with the very weak second ab- 
sorption peak at 4.4-4.5 eV jl9],^0| in polymer thin films. 
Adding weight to this interpretation is the observation 
of a weak but well defined 4.40 eV absorption peak in 
a highly textured film of sexiphenyl (N — 6), orien- 
tated perpendicular to the substrate |l^,|2l|, as well as 
the weak, perpendicularly polarised absorption peak de- 
tected in orientated PFO film in the region 4.2-4.8 eV 
p2[ . In contrast, the l 1 Bi, l state strongly delocalises. 
As can be seen from Fig. pi our DMRG results for the 
^B^ u in the N — 3^. . . , 6 systems practically coincide 
with oligomer data |l7|Jl8[| . For large N the l x B^ u energy 
approaches 3.73 eV in reasonable agreement with the ex- 
perimental peak observed at 3.63-3.68 eV [|l7U2C| 1, or 3.2 
eV 19 1 . Together these results indicate that the first 



(strong) and second (weak) absorption peaks in phenyl- 
based systems are the 1 1 B^ J and l 1 ^^ states, respec- 
tively. 

The third and fourth absorption peaks are polarised 
normal ( 1 i?^ t ) and parallel i}B^ u ) to the long axis, re- 
spectively. In thin films they lie at 5.2-5.3 eV |L9 20 
and 5.7-6.0 eV @|o) respectively. The 1 B~ U state is 
a localised intra-phenyl (Frenkel) excitation which lies at 
6.16 eV in biphenyl ||. To track the position of this state 
as a function of N , we perform calculations targeting a 
number of 1 B^ U states together with the ground state 
l 1 ^.^ and calculate the dipole moments (l 1 A+ -B^) , 
where fi = XiUi (with Xi the x coordinate of the ith 
atom) is the dipole operator along the long molecular 
axis. A typical set of dipole moments, for N = 8, is given 
m Table M. In general, it is found that there are two 1 B lu 
states that are strongly dipole connected to the ground 
state. The first is the l 1 !?-^, i.e., the first absorption 
peak. The second is the localised intra-phenyl excita- 
tion, i.e., the fourth absorption peak. (In the N = 8 case 
this is the ^B^ u state.) The A^-dependence of this high- 
lying, strongly dipole connected 1 B^ U state is plotted in 
Fig.f 

A conspicuous failure of the P-P-P model is its pre- 
diction for the remaining state, namely the lowest lying 
1 B^ U state. The exact solution for biphenyl places this 
state at 6.66 eV, whereas experimentally it is at ca. 5.85 
(and below the 2 1 _B lu state). As shown in Fig. its cal- 
culated energy is 5.9 eV in the long chain limit, 0.6-0.7 
eV higher than the experimental value. Recent work has 
shown that the predicted energy of this state is improved 
if dielectric screening is introduced into the Ohno inter- 
action E3[. 

Fig. |also shows the lowest lying triplet (l 3 ^^) and 
the dipole-forbidden 2 1 state. Electroabsorption stud- 
ies place the 2 1 A+ state at around 4.6 eV [E0U24], around 



0.5 eV below the extrapolated P-P-P result of 5.1 eV. 
(A more recent study jl9| places the 2 1 v4+ substantially 
lower). This discrepancy might possibly be explained in 
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terms of the characteristic red shifts generally observed 
for certain excited states when going from well isolated 
chains to polymers in the solid state. Typical estimates 
for this polarisation or interchain screening shift are ca. 
0.3 eV for the 1 B^ U state and ca. 0.6 eV for the 2 X A+ ||. 
As can be seen from Tabic |y| the 2 1 A+ state has a large 
dipole moment with the l^B^ state, and unlike the case 
for polyenes, it is not pre-dominantly a pair of bound 
magnons, but a particle-hole excitation. (It is usually 
labelled the m}A+ state.) This particle-hole excitation 
is either a 'p'-wave exciton, or the edge of the unbound 
particle-hole continuum. 

To investigate this further, and the nonlinear opti- 
cal properties of poly(phcnylene) in general, we con- 
sider dipole moments between various 1 A£ states and 
the l x Bi u state, as well as between the 2 1 A+ state and 
various 1 B^ U states. The N — 8 values, listed in Ta- 
ble fv|, are representative of the general situation. We 
note that, in addition to the l 1 ^}" and 2 1 A+, another, 
higher lying state, which we denote the k l A+, also has 
an appreciable dipole moment with the l^B-j^. (For the 
N = 8 case k = 5.) We also observed a pattern in the 
(2 1 A+\fi\j 1 Bi u ) values. Namely, the j = 1 state has a 
strong dipole moment with the m}A+, as does the higher 
lying 1 B^ M absorption peak state (localised intra-phenyl 
exciton). In addition, there is another state, lying higher 
still, that has the largest dipole moment with the m 1 . 
We adopt the usual convention of denoting this state as 
the n 1 !?^ (In the N = 8 case n — 7.) 

In Fig. Q we plot a number of 1 A+ and 1 B^ U state tran- 
sition energies as functions of 1/N 2 (^6). We see that 
there are a number of states in the 1 B^ U sector which 
seem to converge to the same energy as the l^-Bj^ in the 
bulk limit. This is the band of 's'-wave charge-transfer 
excitons. Similarly, the m 1 A+ (i.e., the 2 l A+) has a 
number of momentum branches converging to the same 
energy, which is the band of 'p'-wave charge-transfer ex- 
citons. Above this lies the localised intra-phenyl 1 B^ U 
state [p8| . Higher in energy still are the and the 

n B^ u which appear to converge to the same energy in 
the N = oo limit. The strong dipole moments from the 
m 1 A+ to the n l B^ u and the l l B^ u to the k x Aj , and 

the close proximity in energy of the k^A^ and the n x B^ u 
states indicate that these mark the onset of the contin- 
uum of unbound particle-hole excitations. Lying below 
this continuum are the 's'- and 'p'-wave charge-transfer 
excitons and the Frenkel exciton. This hypothesis could 
be checked further by calculating the average particle- 
hole separation for the various states to see if the n^B^ 
is the first unbound state in the 1 B± U sector. The con- 
vergence of the k 1 and the n}B^ u energies to ca. 6.25 
eV as N — > oo would imply a very large binding energy 
(ca. 2.5 eV) for the 1 B± U exciton. However, band states 
are generally expected to be strongly effected by solid 
state screening (a red shift of 1.5 eV has been estimated 
for polyacetylene pj|). Taking such a shift into account 



would bring the n 1 i?^ u and hence the exciton binding 
energy much closer to the results implied by electroab- 
sorption experiments JTs| ] . 

In conclusion, we used a 7r-electron model to investi- 
gate the excited states of poly(phenylene), a model light 
emitting polymer. The transition energies were calcu- 
lated using a new DMRG method. In order to acco- 
modate the large repeat unit in the DMRG calculation, 
a truncated local Hilbert space was constructed for the 
phenyl unit. Using an optimal (or 'natural') basis af- 
fords us convergence in the excitation energies to within 
a small fraction of an eV. For the most part, the calcu- 
lated transition energies agree well with the experimental 
results for biphenyl. By monitoring the evolution of the 
energies with chain length and comparing them to the ex- 
perimental absorption peaks of oligomers and thin film 
polymers, we can assign the four characteristic absorp- 
tion peaks of phenyl-based light emitting polymers. In 
ascending order they are the ^B^, l 1 ^^, l 1 !?^ and a 
high lying, localised intra-phenyl (Frenkel exciton) l 1 i?f tl 
state. A failure of the P-P-P model is the incorrect pre- 
diction that the l 1 ^^ state lies above the Frenkel exci- 
ton |2^] . We have also investigated the finite-size scaling 
of the nonlinear optical states in the longitudinally po- 
larised one- i}B^ u ) and the two- (} A+ ) photon sectors 
and found evidence for two B± u and one 1 j4+ excitonic 
bands lying below the continuum. 
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states in this energy range explicitly in order to clearly see 
higher momentum branches of the localised intra-phenyl 
B^ u exciton. 

[29] Although the positions of these states can be inverted by 
reducing the electron correlation || (the parameters in 
H were determined within CI theory), or by the use of 
more phenomenological, semi-analytic approaches and a 
number of fitting parameters |Bj , our philosophy here has 
been to develop a numerical scheme that can first deter- 
mine the predictions of P-P-P theory (without adjustable 
parameters or uncontrolled approximations), before em- 
barking on the task of attempting to improve it. 



TABLE I. Comparison of DMRG and exact results for the 
TV = 11 phenyl oligomer in the noninteracting (U = 0) limit. 
In the DMRG calculations uia = 280 and tub — 468. 





EQ}A+) 




E(2 l A+) 




DMRG 


-231.0210234 


-228.4828 


-228.331 


-227.195 


EXACT 


-231.0210256 


-228.4832 


-228.328 


-227.213 



TABLE II. Convergence of the lowest two states in the 
(1) for biphenyl (TV = 2). 



sector with m A 



"S 1 


£(lM+) 


E(2 1 A+) 


34 


-63.0095953 


-56.7041530 


92 


-63.0125556 


-56.7137143 


140 


-63.0126802 


-56.7140886 


188 


-63.0126995 


-56.7141368 


260 


-63.0127066 


-56.7141589 


322 


-63.0127078 


-56.7141631 


exact 


-63.0127086 


-56.7141653 
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TABLE III. Convergence of the two-photon (2M+ ) state 
ansition energy with 
for the trimer (N — 3). 



transition energy with the DMRG parameters and trig 



(o) 

m B 


m^> = 34 


m^> = 92 


= 140 


(b) 


420 


5.7626 


5.7498 


5.7486 


574 


5.7578 


5.7449 


5.7437 




692 


5.7562 


5.7433 


5.7421 




952 


5.7544 


5.7415 


5.7405 


(c) 



TABLE IV. Convergence of the 
transition energy with the DMRG 
for the N = 11 system. 



two-photon (2 1 A^) state 
parameters tua and tub 



mB 



34 

92 

140 

202 

268 

346 

418 

468 



140 
204 
322 
420 
574 
692 
848 
952 



34 



85 



rtiA 



130 



niA 



165 



5.38712 
5.23364 
5.21740 
5.20744 
5.20605 
5.20490 
5.20441 



5.35049 
5.17443 



5.17076 
5.16566 
5.16545 
5.16525 



5.36911 
5.19057 
5.16678 
5.15356 
5.15089 
5.14948 
5.14901 



5.36862 
5.18949 
5.16533 
5.15170 
5.14903 
5.14757 
5.14715 
5.14696 



TABLE V. Dipole moments connecting various and 



1 B lu states for the N = 8 system. 



J 


(l^+lAlj 1 ^-) 


O'M+j/ill 1 ^-) 


(2 1 / 


1 


2.85 


2.85 


2.64 


2 


0.68 


2.64 


0.48 


3 


0.19 


0.31 


0.06 


4 


0.11 


0.14 


0.02 


5 


2.52 


1.17 


1.57 


6 


1.03 




1.31 


7 


0.62 




5.06 


8 


0.48 




0.04 




'A '"B A 

FIG. 1. The initial step in developing the optimised 
phenyl block bases used in subsequent DMRG iterations. 
Dashed lines indicate the formation of reduced density ma- 
trices. 



FIG. 2. The steps involved in a DMRG iteration: (a) 
Two system blocks are combined in order to treat systems 
with even N; (b) The middle phenyl unit basis is reoptimised 
in the context of this system size; (c) The system block is 
reoptimised for this superblock; (d) The final energy estimates 
are made, and the system block is augmented to produce the 
system block used in the next DMRG iteration. Dashed lines 
indicate the formation of reduced density matrices. 
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FIG. 3. The transition energies as a function of inverse 
chain length: 1 3 B^ U (pentagons), I 1 B^ u (diamonds), l 1 ^^ 
(stars), 2 1 (squares), l 1 ^^ (circles) and the localised in- 
tra-phenyl (Frenkel) state (triangles). The filled sym- 

bols are the experimental values for biphenyl (see ref. [H]), 
oligomers (N = 3, . . . , 6 [L7 and TV = 6 (jj) and thin 
film polymers (N = oo) [n,[l9]]2i|. The inset shows the 
oligo-phenylene geometry. 
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FIG. 4. The transition energies of a number of A^" and 
1 B^ U states as a function of 1/N 2 , where N is the number of 
repeat units: l 1 B^ u (large, open circles), 2 1 B^[ U (open, down 
triangles), 3 1 B± U (up triangles), 4 1 73^ u (diamonds), the high 
lying localised intra-phenyl 1 B^ U excitation (open squares), 
n l B^ u (pentagons), m = 2 1 (small, solid circles), S 1 ^^ 
(small, solid squares), 4 s - A~^ (x) and k (stars). The large, 
solid circle and solid down triangle show the position of the 
first and second long-axis polarised absoprtion peaks respec- 
tively for sexiphenyl (TV = 6) The dotted lines are to 
guide the eye. 
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